A note on uniqueness in the identification of a spacewise 
dependent source and diffusion coefficient for the heat 

equation 
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We investigate uniqueness in the inverse problem of reconstructing simultaneously a spacewise conduc- 
tivity function and a heat source in the parabolic heat equation from the usual conditions of the direct 
problem and additional information from a supplementary temperature measurement at a given single 
instant of time. In the multi-dimensional case, we use Carleman estimates for parabolic equations to 
obtain a uniqueness result. The given data and the solution domain are sufficiently smooth such that the 
required norms and derivatives of the conductivity, source and solution of the parabolic heat equation exist 
and are continuous throughout the solution domain. These assumptions can be further relaxed using more 
involved estimates and techniques but these lengthy details are not included. Instead, in the special case of 
the one-dimensional heat equation, we give an alternative and rather straightforward proof of uniqueness 
for the inverse problem, based on integral representations of the solution together with density results for 
solutions of the corresponding adjoint problem. In this case, the required regularity conditions on the 
conductivity, source and the solution of the parabolic heat equation are weakened to classes of integrable 
functions. 
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1 Introduction 

Inverse problems of parameter identification in partial differential equations have several important 
applications including thermal prospection of material and bodies, hydraulic prospecting of soil, 
photonic detection of cancer, finding pollution sources, see, for example, [4, 7, 9, 10, 11, 18, 20, 21] 
and references therein. For inverse problems in general it is important to find and specify appropriate 
data such that the set of parameters to be reconstructed are uniquely identifiable. We shall consider 
an inverse problem for the parabolic heat equation, where the additional data is information about 
the solution obtained from a spacewise measurement at the final time. To be more specific, we shall 
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show uniqueness of the identifiability of a pair of functions (a(x), f(x)), representing the spacewise 
thermal conductivity and the heat source, in the parabolic heat equation 



u t - V • (o(i)Vm) = f(x) in Q x (0,T) 

u(x,t) = for (x,t) edttx (0,T) 
0) = h(x) for iGll, 



(1) 



for a given initial temperature h, with the additional temperature measurement g at time t = T, 
i.e. 



We assume that a and / are spacewise dependent real-valued functions and that there exists 
a > such that a(x) > a for all x G f2. This implies that the operator Lu = — V ■ (a(i)V«) 
is elliptic. For the moment we assume that the coefficient, the initial and final conditions are 
sufficiently smooth such that there exist a unique classical solution u(x,t) of the problem (1) and 
that the required compatibility conditions are satisfied. For the precise statement of smoothness and 
other assumptions to guarantee the existence of such a classical solution of the direct problem (1), 
see [15, Theorem 5.2]. Regularity conditions will be further discussed in Section 2 and in the 
Appendix A. 

There are many contributions on uniqueness for the identifiability of coefficients in parabolic 
type equations in the case of lateral overdetermination. Since the literature on this subject is vast 
we suggest the reader to consult [9, 10, 13, 18, 21] and references therein. 

On the other hand, there are only a few papers related to the inverse problem of spacewise 
coefficient identification with given additional measurements at the final time. Indeed, uniqueness 
from final time data for a spacewise dependent heat source was shown in [19] and simultaneous 
determination of a heat source and initial data from spacewise measurements was investigated in 
[12]. Uniqueness from final time data for the identification of diffusion coefficients was shown in a 
recent paper [5], for the one- dimensional heat equation. We point out that the arguments proving 
uniqueness for the one- dimensional setting [5] appear not possible to extend to the multi-dimensional 
case. 

In [21], Carleman type estimates were used to prove uniqueness and stability of a sufficiently 
smooth diffusion coefficient in the heat equation, with measurements given at an intermediate time. 
However, to the authors' knowledge there are no results on uniqueness for the reconstruction of 
both a spacewise dependent conductivity and heat source. Thus, we shall state and prove such 
a uniqueness result building on Carleman estimates in [21]. For the ease of presentation and to 
highlight the usefulness of Carleman estimates for the inverse problem (l)-(2), we shall simplify the 
proof using smoothness assumptions together with some a priori knowledge of the conductivity close 
to the boundary of the solution domain. It is conjectured that these smoothness assumptions can 
be removed but we only indicate possible generalizations rather than give full lengthy and complex 
details. To convince the reader that it is possible to have uniqueness also in spaces of integrable 
functions, we give a rather straightforward proof of uniqueness in the case of a one-dimensional 
solution domain, where the smoothness assumptions are relaxed and more general. This proof is 
based on a recent work [5, Section 5] and involves integral relations obtained from Green's formula 
for the solution together with some denseness properties of the adjoint (heat) equation. This 
technique does not appear possible though to generalize to higher dimensions. 



u(x, T) = g(x) for x G Q , T > . 



(2) 
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We point out that once uniqueness is shown for smooth data, one can use standard approximation 
techniques together with stability results for parabolic equations to get a result for also for non- 
smooth data. 

Note that it is crucial to have only spacewise dependence in the heat conductivity and source 
term; there are examples showing non- uniqueness of the reconstruction of the heat source in the 
case of time-dependent sources, see [2, 6, 8]. A simple counterexample to the uniqueness in the case 
of a time-dependent heat conductivity is presented in [5, Section 5]. 

We point out that an interesting related problem that we do not explore further in this contribu- 
tion is conditional stability estimates for the above inverse problem (l)-(2). In general, for an inverse 
problem, in spite of the ill-posedness conditional stability estimates assure that one can restore the 
stability of the requested (physical) parameters provided they are restricted to some class within 
an a priori bounded set. The conditional stability is not only of theoretical interest but is also of 
importance for the construction of numerically stable solutions. For example, in [3] a stability esti- 
mate gives convergence rates for certain Tikhonov regularized solutions. There are several methods 
in the literature for proving conditional stability, see [9, 21]. The method of Carleman estimates is 
one possibility for obtaining conditional stability. Thus, it is possible to obtain such stability based 
on the results presented in the present paper. We do not go into the details here, it is deferred to 
future work. The reader can further consult [9, 18, 21] and the references therein. 

For the outline of this paper, we show a uniqueness result for the above inverse problem using 
Carleman estimates, see Section 2 and Theorem 4. In Section 3, we address global uniqueness for 
the identification of the heat source and conductivity in the one-dimensional case, using density 
arguments for solutions of the corresponding adjoint problem. The arguments in Section 3 weakens 
the smoothness assumptions of the spacewise source and coefficients in the inverse problem (1)- 
(2). In the final section, we draw some conclusions and discuss some possible generalization of the 
presented uniqueness results. For the sake of completeness, in the Appendix A, we prove that for a 
given class of initial data, the solution of the heat equation at the final time satisfies the required 
assumptions stated in Section 2. 

Notation: We finish the introduction stating some notation that we use: The set Q is an open 
and bounded subset of M n , with the boundary dQ being at least Lipschitz smooth. By L P (Q) for 
1 < p < oo, we denote the usual space of p-integrable functions on Q with the usual norm IHI^™- 

The space L°°(Q) is the standard L°°-space. We denote by W k ' p (Q) the standard Sobolev space on 
Q with generalized derivatives of order < k in L P (Q). In particular, for p = 2 we have the Hilbert 
spaces H k (Q). Moreover, since dQ is Lipschitz, the trace of a function in H 1 ^) to the boundary is 
well-defined. 

Let T > be fixed and define the measurable function u(-,t) : (0,T) — > X, where X is a 
Banach space. We denote by C([0, T];X) the space of continuous mappings u(-,t) with the usual 
norm and by L 2 ((0, T); X) the space of functions such that 



Jo 

We also assume that a and / are sufficiently regular spacewise real-valued functions with < 
a < a(x) for every igSI. Therefore, the differential operator 

Lu = V(o(i)Vu) 
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is elliptic. Since, the operator — L generate a contraction semigroup, there exists a unique solution 
u of (1) with 



2 Local uniqueness for the spacewise source and heat conduc- 
tivity in (l)-(2) 



As mentioned, we do not strive to obtain the most general result, but we only wish to convince 
the reader that there can be at most one spacewise dependent heat conductivity coefficient and 
one spacewise dependent heat source that together satisfy the given final time data, i.e. that solve 
the inverse problem (l)-(2). Therefore, we assume that the given data and the solution domain 
are sufficiently smooth such that the required norms and derivatives of the coefficients, sources 
and solution of (l)-(2) exist and are continuous throughout the solution domain. For the precise 
statement of smoothness and compatibility conditions for the parabolic heat equation, see [15, 
Theorem 5.2]. Once uniqueness is shown for smooth data, one can use standard approximation 
techniques together with stability results for parabolic equations to get a result for non-smooth 
data. Note that there are now many solvability results and estimates for parabolic equations with 
very general coefficients, see further [15, 14]. 

For our proof of uniqueness we use local Carleman estimates and follow [21, Section 6]. To 
simplify the presentation further and to avoid cut-off functions and global Carleman estimates, we 
assume that the diffusion coefficient is known in a region Q — D, where Dccfl, and < \D\ < oo. 
In other words, 

Assumption 1. The diffusion coefficient a(x) in (1) is known for every x G Q \ D, where D C Q 
with dD sufficiently smooth and d(D, dfl) > 7 > 0. 

This is a reasonable assumption in applications, since the material (body) Q might be coated or 
layered and the physical properties of the outer layer is known. The above assumption forces any 
two solutions of the inverse problem (l)-(2) to be equal in fi \ D, i.e. the difference has compact 
support in Q. More precisely, we have the following lemma: 

Lemma 2. Let u = u(a, fx) and v = v(b, f 2 ) be solutions of (l)-(2), with spacewise heat conduc- 
tivities a and b and spacewise heat sources j\ and $2, respectively. Assume that a = b in Q — D. 
Then the difference w = u — v is such that w = in Q\D. Moreover fi = f 2 in Q — D. 

Proof. By linearity of (1), w satisfies 
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dt < 00 . 



wt-V ■ (a(x)Vw) = V ■ ((a - b)Vv) + {f x - j 
w{x,t) = for (x,t) edttx (0,T) 
w(x, 0) = w(x, T) = for x G fl 



f 2 )(x) in Ox(0,T) 



(3) 
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From the assumption that a = b in Q — D we have, 

w t -V- (a(x)Vw) = (fx - f 2 )(x) in Q - D x (0, T) 

w(x, t) = for (x, t) G dtt x (0, T) (4) 
w(a;,t) = ^ (a;) for (x,t) e dD x (0,T) 
iy(x, 0) = w(x, T) = for x G f2. 

Therefore, using standard parabolic theory we conclude that there exist only one w G C([0, T], Hq(Q) x 
H 2 {Sl) solution of (4). 

Let we define the following problems 



w 



(1) - V • (a(x)Vw (1) ) = U'i ~ h){x) in Q - D x (0,T) 



t» (1) (x, t) = for (i,i)e90x (0, T) (5) 
io (1) (a;,t) = for e dD x (0,T) 
w (1) (x,0) = w; (1) (x,T) = for iGd, 

Since the source is in L 2 (f2), we have from [19, Theorem 2] that there exist a unique solution it/ 1 ) 
of (5). 

Therefore, taking the difference it/ 2 ) = if — it/ 1 ), we find that it/ 2 - 1 exists and satisfies 



w 



(2) - V • (a(x)Vn; (2) ) = in Q - D x (0, T) 



«/ 2 ) ( Xj t) = for (x, t) G dQ x (0, T) (6) 
w (2 \x,t) = g (x,t) for e dD x (0,T) 
w {2) (x,0) = w (2) (x,T) = for iGfl. 

We start by considering the solution it/ 1 ) of (5). From [19, Theorem 1] given a final condition 
at time T > there is at most one spacewise dependent source giving rise to this final time value 
for the heat equation with homogeneous boundary and initial condition. Thus, since it/ 1 ) is zero at 
t = T and has homogeneous boundary and initial condition we conclude that f\ — f% = 0, this in 
turn implies that it/ 1 ) = in f2 — D x (0, T). 

We then consider the solution it/ 2 ) of (6). If the boundary data go is identically zero it is clear 
that it/ 2 ) = 0. Therefore, assume that go is not identically zero. Without loss of generality, we 
can assume that g (x,t) > at a point (x,t) of the boundary dD x (0,T) and since go is a ^ least 
continuous it follows that g (x,t) > in E x (ti,t 2 ), with E being a surface patch of dD. This in 
turn implies that it/ 2 )(a;,t) > for (x,t) sufficiently close to this open set. Applying [7, Theorem 
9.2] it follows that it/ 2 )(x,s) > for every t < s < T contradicting the given final condition 



w 



(2), 



x,T) = 0. Therefore, go = implying that it/ 2 ) = 0. 
Thus, since w = u/^+u/ 2 ) we have w = and therefore the assertion of the lemma is proved. □ 



Note that, if the measurements are given at an intermediate time < to < T, the conclusions 
of Lemma 2 may not be true. 
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2.1 Some additional assumptions and admissible solutions 



In order to proceed to the next step in the proof of uniqueness of the spacewise dependent pair 
(a,/) satisfying (l)-(2) we need some additional assumptions. We assume that the admissible set 
of unknown elements is 

A:={(aJ)EC 2 (U) : a(x) >a>0,i€fl, \H c ^(n) + WfW&M ^ M > > ( 7 ) 

where M > is an arbitrary fixed constant. Moreover, we assume that g and h are smooth enough 
such that we can take the t-derivatives of w = u — v for the time points needed, see [15] for the 
details on the required smoothness assumptions on g and h for this to be the case. Note that in 
Section 3 and in Appendix A we will provide more details on the regularity of the parameters and 
the input data. Moreover, we also assume that dQ is smooth enough such that dQ C {x± = 0}. 
In fact, the boundary can be covered by a finite number of surface patches mapping to this region 
and it is therefore enough to consider the estimates in this half-space. The general results can be 
obtained in standard way for partial differential equations by using a partition of unity argument. 

We remark that since the coefficients and the source term are time-independent, the solution of 
(1) is analytic in time [9, 20]. Moreover, the compact support of w in Q guaranteed by Lemma 2 
implies the follows boundary conditions 

(x,t) Gfflx (0,T), 

(x, t) G dQ x (0, T) , (8) 
(x,t) Gfflx (0,T), 

where t] is the outer normal vector and r is the tangential normal vector at dQ. Moreover, 
-7^w t (x,t) = 7](x) ■ Vu>t(x, •) and -^w t (x,t) = t(x) ■ Vw t (x, •) are the normal and tangential deriva- 
tives of w t for x G dQ. 

We shall present a uniqueness result for the inverse problem (l)-(2) based on local Carleman 
estimates. A crucial step for such estimates is the construction of suitable weight functions, <p(x,t) 
and (3(x, t). For a given domain D we can apply the arguments in [21, Section 5.1] or the arguments 
in [18, Section 2], with a special choice of a function d(x) such that 

<p(x,t) = e x ^ x ' t) (9) 

and /3(x,t) = d(x) + e(t), where max tg [ 0i T] f(x,t) = <p(x,T). Moreover, we can construct <p and (3 
with 

Q := {(x,t) : Xi > , (p(x, t) > e~ xs } = {(x, t) : X\ > , fi(x, t) > -5} , (10) 

and Q H {t = T} = D. In (9) and (10), A, 5 > are some fixed constants being sufficiently large. 
Let we give a very simple example for choice d(x) and e(t) form [21, Section 6]. 

Example 3. Set x' = (x2,-,x n ) and x = (xi,x') G IR n . Let we assume that D = D(8) : = 
{(x, x');0 < x\ < — \x'\ 2 / , ~f + 5/7} and moreover D(45) C Q, for some 7 > and 5 > 0. Define 
d(x) := — 7x1 — \x'\ 2 ande(t) := —(t—T) 2 . Therefore, we have by (9) that max tG [o,T] ^{x^t) = ip(x,T) 
and that Q = Q(5) satisfies Q D {t — T} = D. 



w t {x,t) = 

d 

—w t (x,t) = 
—w t {x,t) = 
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As mentioned earlier, since we have spacewise dependent coefficients, the solution of (1) is 
analytic in time. Therefore, the solution can be extended beyond the final time T, see [20, Section 
3] and T can therefore be considered as an interior point, as is required in the derivation of the 
results in [21, Section 6]. 

We remark that one can try to identify the most general set of conductivity coefficients having 
minimal regularity assumptions together with minimal regularity of the boundary of Q for which 
the weight functions <p and (3 do exist. However, as pointed out earlier this is not the main aim of 
this study, we shall only present a proof in the case of smooth and regular solutions and domains. 
For a discussion about more general function spaces and domains for which this derivation can hold 
true, see [18, 21] and references therein. 

We also remark that since w has compact support (see Lemma 2), we do not need to use cut-off 
functions. Therefore, many of the calculations and terms in [21, Theorem 6.1] can be dropped 
and the steps in the proof become easier to follow. However, since we have two unknowns, some 
challenges and adjustments do remain, and we present the steps below. 



2.2 A uniqueness proof 

Let u = u(a, fx) and v = v(b, f 2 ) be solutions of (l)-(2), with spacewise heat diffusions a and b 
and spacewise heat sources fi and f'2, respectively, and put w = u — v. Denote by z := w t ; a well- 
defined quantity due to the smoothness assumptions above and since we only work with classical 
solutions. Given the analyticity of w in time and the assumption that the unknown pair (a, /) is 
time independent, it follows that z satisfies 

z t -V ■ (a(x)Vz) = V ■ ((a-b)Vv t ) (11) 

with homogeneous boundary, initial and final conditions. Then, [21, Theorem 3.2] implies 

/ \-J2 \did jZ \ 2 + s\\/z\ 2 + s 2 \z\ 2 \e 2sip dxdt<C [ |V • ((a - b)Vv t )\ 2 e 2s<f dxdt (12) 
J Q\ S i, j= i J JQ 

for all sufficiently large s > 0. Note that the integral over the boundary in [21, Theorem 3.2] is 
identically zero due to (8). 

Now we shall derive an estimate of the right-hand side of (12). A direct calculation of V ■ ((a — 
b)Vv t ) implies that 



|V • ((a - b)Vv t )\ 2 < (|V(a - b)\ 2 + \a - b\ 2 )(\\/v t \ 2 + |A^| 2 ) . (13) 



Hence, 



/ \V-((a-b)Vv t )\ 2 e 2 ^dxdt<C(\\Vv t \\l + \\Av t \\ 2 J [ (|V(a - b)\ 2 + \a - b\ 2 )e 2s * dxdt . (14) 

JQ JQ 

Since u, v are assumed smooth enough (see [15] for necessary requirement on the data) we have that 
(\\VvtWl + \\Avt\H) < 00. Therefore, C = (||V^||^ + ||A^||^) is finite. 
Combining (12) and (14) imply 

7 (a-b)\ 2 + \a-b\ 2 )e 2s * dxdt (15) 



/ ( J2 \did 3 z\ 2 + s 2 \Vz\ 2 + s 3 \z\ 2 ) e 2sLp dxdt<C [ s(|V( 
J Q Vu=i / J Q 
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for all sufficiently large s > 0. 

Given the final time measurement (2), equation (3) gives, at t = T, 

V ■ ((a - b)Vg(x)) = w t (x, T) + (fx - f 2 )(x) (16) 
V(V • ((a - b)Vg(x))) = Vw t (x, T) + V(A - f 2 )(x) . 

It follows from a straightforward manipulation of (16) together with integration that 
/ (|/i-/ 2 | 2 + |V(A-/ 2 )| 2 )e 2 ^' T) ^ 

JD 

+ [ (\V-((a-b))Vg(x)\ 2 + \V(V-((a-b)Vg(x)))\ 2 )e 2s ^ T Ux 

JD 

< [ (\w t (x,T)\ 2 + \Vw t (x,T)\ 2 )e 2s ^ T Ux (17) 
jd 

+ 2 / \(V-((a-b)Vg(x))(f 2 -f 1 )(x))\e 2s ^ T Ux 
Jd 

+ 2 / \(V(V-((a-b)Vg(x)))V(f 2 -h)(x))\e 2s ^ T Ux. 
Jd 

We put 

A = 2 [ \(V-((a-b)Vg(x))(f 2 -f 1 )(x))\e 2s ^ T Ux 



D 



and 



B = 2 [ \(V(V-((a-b)Vg(x)))V(f 2 -f 1 )(x))\e 2s ^ T Ux. 
Jd 



Note that 



\\djW t e s % 1{Q) < ( {^^O^ + ^w^ + s^wAe^dxdt. (18) 

jQ \i,j=l J 

Applying the trace theorem [1] in Q and noting that Q D {t = T} = D we have 

/ (K(*,T)| 2 + |V^(x,T)| 2 )e 2s ^ r ) dx < Cj2\\d jWt e s m i{Q) . (19) 
Jd ._ n 



3=0 

Now, (15), (17), (18) and (19) yield 



/ (|/i-/2| 2 + |V(A-/ 2 )| 2 )e 2 ^ T )rfx 
Jd 

+ f {\V-((a-b))Vg(x)\ 2 + \V(V-((a-b)Vg(x)))\ 2 )e 2s ^ x > T Ux 

JD 

<C f s(|V(a -b)\ 2 + \a- &| 2 ) e 2s ^ dxdt + (A + B). (20) 
JQ 

The next step in the proof is an estimate of the second integral in the left hand-side of equa- 
tion (20). 



Lemma 4. [21, Lemma 6.1] Let (a — b) E H 2 (tt) such that \a — b\ = |V(a — fe)| = ondVt. Assume 
that 

n 

1 d 1 g(x)+2j2(d j g)(x)<0 xeQ, (21) 
i=2 

dig(x)>0, xedQ. 

Then there exists a constant C > such that 
! (|V(a-6)| 2 + |a-6| 2 )e 2s ^ T) dx 

<C [ (\V ■ ((a - b))Vg{x)\ 2 + \V(V ■ ((a - b)Vg(x)))\ 2 )e 2s ^ x ' T) dx (22) 
Jn 

for all sufficiently large s > 0. 

Note that, since we assume that a = b in Q — D, the estimate (22) holds over D. This in turn 
using equation (20) and Lemma 4 give 

(l/i - / 2 | 2 + |V(/i - f 2 )\ 2 )e 2s ^ T) dx + / s 2 (|V(a - b)\ 2 + \a - b\ 2 )e 2s ^ x ^ dx 

D JD 

< C [ s(|V(a - 6) | 2 + \a - b\ 2 )e 2s ^ t] dxdt +{A + B) (23) 

JQ 

for all sufficiently large s > 0. Since tp(x, t) < tp(x, T) , x G (fi) , < t < T, we have 

(l/i - / 2 | 2 + |V(A - / 2 )| 2 )e 2s ^ T > rfx + / s 2 (|V(a - b)\ 2 + |a - 6| 2 )e 2s ^ T > dx 

< CT f s(|V(a - 6) | 2 + |a - 6| 2 )e 2s ^' T) dxdt + (A + B) (24) 

for all sufficiently large s > 0. 

The final step in the proof of uniqueness is to obtain an estimate of the quantity A + B in (24) 
in terms of the coefficients of the first equation in (3). We assume that 

IIVCA^*))^ + UAyOiOL + HVylL < oo . (25) 

Using the Young inequality with e > we have 

A = 2 f \(V-((a-b)Vg(x))(f 2 -f\)(x))\e 2s ^ T Ux<- f |(V ■ ((a - b)V g{x))\ 2 e 2s ^ T) dx 

JD £ JD 



+ e / \{f 2 -h){x))\ 2 e 2s ^ T Ux 

JD 

and 

B = 2 [ |(V(V-((a-6)V^(x)))V(/ 2 -A)(x))|e 2 ^ T )dx<- / |(V(V ■ ((a - b)Vg(x)))\ 2 e 2s ^> T Ux 

JD £ JD 

+ e [ \V{f2-h){x))\ 2 e 2s ^ T Ux. 



D 



Using the similar estimate as before for |(V ■ ((a — b)Vg(x))\, see (14), in combination with 
assumption (25) we have 



A<- (\V{a-b)\ 2 +\a-b\ 2 )e 2s ^ T) dx + e / \(f 2 - fi)(x))\ 2 e 2s ^ T) dx. 
£ Jd Jd 

Similarly, a direct calculation of (V(V • ((a — b)Vg(x))) using the assumption (25) imply 

/ |(V(V • ((a - b)Vg(x)))\ 2 e 2s ^ x ' T) dx < C [ (|A(a - b)\ 2 + |V(a - b)\ 2 + \a - b\ 2 )e 2s ^ x ' T) dx. 
Jd Jd 

We then choose < e < 1/2. Using the above estimates for A and B in (24) and put together 
this give 

/ (|/i-/2r + |V(/ 1 -/ 2 )| 2 )e 2 ^ T )^+ / s 2 (\V(a-b)\ 2 + \a-b\ 2 )e 2s ^ T Ux 
Jd Jd 

<CT [ [(a + l)(|V(a - b)\ 2 + \a - b\ 2 ) + |A(a - b)\ 2 }e 2s ^ x ' T) dxdt (26) 
Jd 

for all sufficiently large s > 0. 

For s > sufficiently large the terms in the right-hand side can be absorbed by the corresponding 
terms in left hand-side. Indeed, by assumption, |A(a — b)\ 2 < M 2 and therefore it is enough to 
choose s > 0, such that s 2 (|V(a - b)\ 2 + \a - b\ 2 ) > CT(s + l)(|V(a - b)\ 2 + |a - 6| 2 ) + M 2 . Thus, 
with such a choice of s > 0, we conclude that the second integral in the left-hand side is identically 
zero, i.e. a = b also in D. The inequality (26) then implies that also the sources are equal, i.e. 
fi — f2- Since we concluded that the coefficients are equal one can alternatively use [19, Theorem 
1] to obtain uniqueness of the sources. 

The obtained results can be summarised as follows. 

Theorem 5. Let the spacewise conductivity coefficient and heat source (a, f) £ A, where A is given 
by (7) , with the coefficient a satisfying Assumption 1. Moreover, assume that the initial condition 
h and the final time measurement g are regular enough such that the corresponding solution u(x,t) 
of (1) satisfies 

Moreover, assume that the final time condition g satisfies the Assumption (21) in Lemma 4 and 
(25). Then the inverse problem of identifying {a(x), f(x),u(a, /)} in the heat equation (1) for a 
given additional measurement g(x) = u(x,T) has a unique solution. 

In the Appendix we shall prove that there exist conductivities and sources which can generate 
a final time value satisfying (21), i.e. the inverse problem (l)-(2) with (21) imposed will have a 
solution for some data (and this solution is unique according to the above theorem). 

We could now dwell into lengthy calculations on the uniqueness of the inverse problem (1)^(2) 
under less regularity assumptions on the coefficient and data. There are indeed many generalizations 
of Carleman estimates for parabolic equations that one could potentially use in this case similar to 
the arguments in [18, 17, 20, 21]. However, we prefer not to enter into these technicalities. 

Let us though briefly discuss some alternative estimates in the above proof of uniqueness that 
could potentially reduce the imposed smootheness assumptions. We used the estimate (13) to obtain 
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an upper bound of the right-hand side of (12). However, there are alternative estimates that can 
be used. For example, in [18, Subsection 3.2] the right-hand side of (12) is estimated by 



Note that, with the above estimate, equation (15) follows with less assumptions on the regularity 
of the coefficients as well as of (1). Moreover, in [18, Theorem 3.1] a Carleman estimate (for L°°(Q) 
coefficients) relating the L 2 -norm of the coefficients with the right-hand side of equation (16) is 
given. In other words, [18, Thoerem 3.1] implies that one do not need to differentiate as high 
in time as we did in the second equation of (16). These two factors together would, most likely, 
improve and weaken the smoothness assumptions in Theorem 5. 

We point out that the Assumption 1 was imposed merely for technical reasons to simplify the 
presentation of the Carleman estimates and to avoid cut-off functions. To convince the reader that 
this assumption can be removed and to further highlight that our assumptions are far from the most 
general ones, we discuss in the next section the uniqueness for a more general class of spacewise 
parameters in the one-dimensional parabolic heat equation using a completely different technique, 
which does not lend itself to higher dimensions though. 

3 Uniqueness of the spacewise heat conductivity and heat 
source in a one-dimensional heat equation 

In this section we show the uniqueness under weaker smoothness assumptions compared with 
Section 2 of the identification of the spacewise pair of coefficient and source (a,f) in (1), with 
additional measurements given by (2). In particular, the Assumption 1 will not be used. However, 
we only consider the case with Q = (0, L) C K., i.e. the spacewise solution domain in (l)-(2) is 
one-dimensional in space. 

The derivation of the uniqueness result is based on a completely different technique than 
Carleman estimates. Indeed, the technique is based on results that relate the uniqueness of the 
inverse identification problem to the density in certain function spaces of solutions of the corre- 
sponding adjoint problem. With this approach the assumptions on the smoothness of the unknown 
pair (a, /) is reduced to a more general class. Moreover, assumptions on the smoothness of the 
input and measured data are determined only by extracting a certain differential dependence and 
that there exists a solution for the corresponding direct problem. 

We start by assuming the following regularity conditions for the parameter, the source, the initial 
condition and the measured data in the inverse problems (l)-(2); compare with the assumptions in 
Section 2. 

Assumption 6. We assume that the heat conductivity a £ L°°((0,L)) and that there exists a > 
such that a(x) > a for all x £ (0, L) . Moreover, it is assumed that a(-) is continuous in [0, e) and in 
(L — e, L] for any fixed e > 0, the heat source f £ L 2 ([0, L]), the initial temperature h £ -ff 1 ([0, L\), 
and that the additional final time temperature measurement g £ -ff 1 ([0, L]) . 

Since h,g £ iJ 1 ([0,L]) we further need to assume that the following matching conditions 




Q 
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(compatibility) are satisfied in x = and x = L 



-{a{0)h x (0)) x = f(0), 

-(a(L)h x (L)) x = f(L), (27) 
-(a(0)g x (0)) x = /(0) , 
-(a(L)g x (L)) x = f(L). 

For the existence and regularity of a solution for the corresponding direct problem with the pair 
(a(x), f(x)) satisfying the conditions on regularity stated in Assumption 6, we have: 

Lemma 7. Let the Assumption 6 and the matching conditions (27) hold. Then, there exists a 
unique solution u(x,t) G H 2,1 ([0,L] x [0,T]) of the boundary value problem (1). Moreover, there 
exists a constant M such that \\ u \\ C ([o l]x[ot]) — 

Proof. The existence and uniqueness follows directly from classical results on parabolic partial 
differential equations, see for example [16, 15]. Now, from the Sobolev embedding Theorem [15, 16, 
1] we have that u(x, t) G C([0,L] x [0,T]). The uniform boundedness follows from the maximum 
principle for parabolic equations [15] together with the assumed smoothness of the boundary, initial 
and final data. □ 



3.1 Uniqueness of a solution of the inverse problem: 1-d case 

The steps for proving uniqueness of the identification of the pair of parameters {a(x),f(x)} for 
given initial and final data in (l)-(2) are outlined below: 

Assume that u = u(a, f\) and v = u(b, f 2 ) are two solutions of (1) with additional data (2). As 
in the previous section let w — u — v. Then, w satisfies 

w t - (a(x)w x ) x = ([a(x) - b(x)]v x ) x + (A(x) - f 2 {x)) in (0, L) x (0, T) (28) 

with homogeneous initial, boundary and final conditions. 

For the proof uniqueness in the one-dimensional case we shall invoke the adjoint problem of 
(28), that reads as 

A + (flWi), = in (0,L) x (0,T) 
V>(0, t) = ip(L, t) = for t E (0, T) (29) 
ip(x, 0) = for x e (0, L) 
tp(x,T) = ji(x) for x e (0, L), 

where /i(x) is an arbitrary function in Cq [0, L]. 

For properties of solutions of the adjoint equation (29) we have: 

Lemma 8. Let the Assumption 6 hold. 

i) For any function n{x) G Cq[0,L], there exists a unique solution i/;(x,t;v) G C 1 ((0, T); C 2 (0, L))fl 
C([0,L] x [0,T]) of (29). 
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ii) For any function n(x) G Cg [0, L], the following relation holds 

! j tf>(x,t;[i(x))F(x,t)dxdt = 0, (30) 
Jo Jo 

where w is a solution to (28) with right-hand side F(x,t). 

Hi) For fi(x) ranging over the space Cq[0,L], the corresponding range of i/)(x,t; fi(x))\ t=T is every- 
where dense in L 2 [0, L] for any time t = r, < r < T. 

iv) Given that 



i/j(x, t; /i(x))$(x, t) dxdt = 
for n(x) ranging over the space Cq[0,L], then 

$(x,T) = 0. 



o Jo 

2\ 



Proof. Item i) is a well-known result for parabolic equations, see, for example, [16, 15]. 
Item ii) follows immediately by multiplication of (32) by ip and integration by parts. 
Item iii) and Item iv) are consequences of [5, Lemma 2-3, pg 318]. 

□ 

We now have the required results in order to prove the main step in the uniqueness argument for 
the conductivity function a(x) by adding some additional assumptions on the final data in (l)-(2). 

Theorem 9. Let the Assumption 6 and the matching condition (27) hold. Moreover, assume that 
\g x {x)\ > for all < x < L. Then the inverse problem (l)-(2) has a unique solution {u, a, /} with 
the conductivity a G L°°(0,L) ; the heat source / G L 2 (0,L), and temperature u, with \\u\\ l < oo. 

Proof. Since we have time-independent coefficients and source, it follows that the solution v to 
(1) with additional data (2) has derivatives of all orders with respect to t [9, 20] and this in 
turn implies that w has derivatives of all orders with respect to t. Moreover, due to the Sobolev 
imbedding theorem, we can assume that w in (28) is at least continuous; for simplicity we assume 
that pointwise evaluation makes sense for the coefficient and source. 

Note that at x = and at x = L using the assumptions and matching condition (27) imply that 

a(x) = b(x) for x = and x = L . (31) 

In order to prove that a(x) = b(x) for < x < L we apply Lemma 8 Item ii) in combination 
with iv) in (30) to get 

[{a{x) - b{x)){v(x, T)) x ] x + {h{x) + f 2 {x)) = 

for < x < L. Using this in the first equation in (28) in combination with w(x, T) = for every x 
in [0,L], we conclude that w t (x,T) = 0. 
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Define z = Wt- Similar to the previous section we have that z satisfies 

z t - (a(x)z x ) x = ([a(x) - b(x)](v x ) t ) x in (0, L) x (0,T) 
z(0, t) = z(L, t) = for t e (0, T) (32) 
z(x, 0) = e^x) for x e (0,L) 
z(x,T) = for x e (0,L). 

Splitting this problem into two, one with zero right-hand side and with initial condition 9\ and one 
with the given right-hand side and zero initial condition, following the proof of Lemma 2 one can 
conclude that the solution to the first one is identically zero, i.e. 9\{x) = 0. Therefore, since z 
satisfies a problem of the same kind as w we can again apply Lemma 8 Item ii) in combination with 
iv) in (30) to conclude that 

[(a(x)-b(x))(v t (x,T)) x ] x = 0. 

Using this in the first equation in (32) in combination with z(x,T) = for every x in [0,L], we 
conclude that z t (x,T) = 0, i.e. w tt (x,T) = 0. Continuing this, putting z 1 = z t and deriving the 
problem for z\ and applying the similar reasoning, i.e. Lemma 8 Item ii) in combination with iv), 
we find that Wut{x,T) = 0. Further continuing this it is possible to prove that (d ( t k) w)(x,T) = 
for A; = 0,1,2,... and the same holds at t — 0. From this and strong unique continuation results 
for parabolic equations, we conclude that w(x,t) = for [0, L\ x [0,T]. This in particular implies 
that z = in [0, L] x [0, T}. From the first equation in (32) we then have 

{(a(x) - b(x))(v t (x,t)) x ] x = 

for every (x, t) in [0,L] x [0,T]. Integrating first with respect to x using that a(0) = 6(0), we find 
that 

(a(x) - b(x))(v t (x,t)) x = 0. 
Since the coefficients are independent of time, we write this as 

[(a(x)-b(x))(v x (x,t))] t = 0. 

Integrating with respect to time using a(0) = 6(0), and then putting t = T we obtain 

(a(x) - b(x))g x (x) = 

for < x < L. From the assumptions on g we can conclude that a(x) = b(x) also for < x < L. 

The final step in the uniqueness argument is the proof of unique identifiability of f(x) in (l)-(2). 
Since w = and a = b, we have from (28) that 

/i(z)-/ 2 (x) = 

for every x in [0, L], i.e. f\ = / 2 and the theorem is proved. □ 

We finish remarking that the main uniqueness argument is related to the thermal diffusivity 
coefficient a. However, the arguments that we have presented in this section appear not possible to 
generalize to dimensions in space higher than one. This is due to the fact that we would obtain an 
equation where the divergence of an element is zero. However, we can not conclude that the given 
element would then be a constant, as the simple example (x, —y) shows. 
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4 Conclusions and possible generalizations 



In this paper, we proved uniqueness for the inverse problem of simultaneously identifying a spacewise 
heat conductivity and heat source for a given final time measurement. The main result is based 
on local Carleman estimates for parabolic problems following [21]. We did not strive for the most 
general result, but only aimed at convincing the reader that there can be at most one coefficient 
and source that satisfy a given final time condition. Therefore, we assumed rather strong regularity 
on the unknown parameters and worked with classical solutions. However, these can be relaxed 
and it was indicated how to adjust the proof. For conductivities and sources in spaces of integrable 
functions an alternative proof of uniqueness in the one-dimensional case was given, where many 
of the assumptions on smoothness were weakend. Unfortunately, these arguments for the one- 
dimensional case can not be generalized to dimensions higher than one. The one- dimensional 
case does motivate the uniqueness of the inverse problem under less regularity assumptions on the 
parameter spaces; to present such a result in higher dimension is deferred to future work. 

We remark that there are actually many strong results on Carleman estimates for parabolic 
equations with very general assumptions on the smoothness of the parameters [18, 17]. Moreover, we 
used local Carleman estimates in order to simplify the proof of uniqueness in the higher- dimensional 
case. However, we may use global Carleman estimates to avoid the assumption that the heat 
conductivity is known close to the boundary of the body Q, see [13, 17, 18, 20, 21] for a general 
overview on the subject. 

Therefore, the authors' conjecture is that there exists uniqueness of the spacewise heat conduc- 
tivity and heat source for a given final time additional measurement, under weaker assumptions on 
the parameter space. It will be addressed in a future works together with a regularization method 
for the inverse problem. 
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A Appendix 

In this appendix, we show that the assumption (21) on the final time measurement (2) required in 
Section 2, can be satisfied, i.e. the inverse problem (l)-(2) with (21) imposed can have a solution 
for a given initial data being sufficiently regular and with a localized heat source. 

The first result is on the regularity of the solution of (1). For simplicity, we assume that all 
compatibility conditions are satisfied and that Q is an open, bounded and regular subset of M. n . We 
do not present a proof of the following theorem, since it follows from standard regularity estimates 
for parabolic equations. Moreover, since the coefficient and the source are spacewise dependent, the 
solution is analytic in time. For details see, for example, [16, 15]. 

Theorem A.l. Let the coefficient and the source belong to the admissible set A. Then, for any 
h E H k (tt), there exists a unique solution u E C([0, T]; H k+1 (£l)) D C((0, T); Hq(Q) D H k+1 (il)) n 
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C\(0,T);H k+1 {Q)) of the parabolic equation (1). 

Note that the regularity above is far from optimal. 

Lemma A. 2. Let k > max{3,n/2}. Then the assumption about the finiteness of the constant C 
in (15) and assumption (25) are satisfied for the solution of (1) with additional data (2). 

Proof. This follows from Theorem A.l and the continuous embedding of H k (Q) (1 Hq(Q) in C(Q), 
[!]■ ' ' □ 

Now, we verify Assumption 21 in Lemma 4. We remark that Lemma 4 is a corollary of [21, 
Lemma 6.2] that reads as follows 

Proposition A. 3. Let the solution domain fl be as above. Put ( = a — b. Consider the first-order 
partial differential equation 

(Po()(x) = VC(x) ■ Vg(x) + C(x)Ag(x) 

where g G H k (Q), for k as in Lemma A.2. 

If the Carleman weight function tp(x,T) G C 1 ^) satisfies 

Vg(x) ■ V<p(x,T) > 0, xgH, (A.l) 

then the Carleman estimate (22) is satisfied. 

Therefore, we only need to prove that there exists a Carleman weight function ip, a set of initial 
data and a heat source such that (A.l) is satisfied. However, by construction of the Carleman 
weight function, we have 

V<p(x,T) = Vd(x)e x ^ x ' T) . 
Hence, is enough to guarantee that there exist d(x) such that 

Vg(x) ■ Vd{x) > , xeQ. (A.2) 

It is verified in the following lemma. 

Lemma A. 4. Let s > and let O be any open set of Q. There exists d G C{Q) with d\gn = and 
\Vd(x)\ > e for x G fl, such that for any h G L 2 (Q) there is a sufficiently smooth source term f , 
having support in O, with 

\\\7u(f,h)(T)-\7d\\ Lx{n) <e/V2, (A.3) 

where u(f,h) is the solution of (1) with source f and initial data h. Moreover, (A.2) holds. 

Proof. Since we do not use global Carleman estimates, we can consider d to be zero on dQ. For the 
existence of such a function d, see for example [18, Section 2]. The density argument then follows 
from [22, Proposition 1.1], see also [18, Corollary 3.1]. 

Hence, from g{x) = u(f, h)(T), the regularity of g(x) and d(x) and the estimate (A.3), we have 

e 2 /2 > (esssup \Vg(x) - Vd(x)\) 2 > \Vg(x) - Vd(x)\ 2 = |V#(x)| 2 - 2Vg(x) ■ Vd(x) + \Vd(x)\ 2 . 

Therefore, since \Vd\ > e it follows that 

2Vg(x) ■ Vd(x) > \Vg(x)\ 2 + \Vd(x)\ 2 - e 2 /2 > \Vd(x)\ 2 - e 2 /2 > e 2 /2 . 

□ 



16 



The condition (A. 2) can then be stated as (21) via a suitable transformation, see the proof of 
[21, Lemma 6.1]. Thus, we have shown that there are conductivities and sources that can generate 
a final time value such that (21) holds, i.e. the inverse problem (l)-(2) with (21) imposed will have 
a solution for some data (and this solution is unique). Again, we have not investigated optimal 
conditions and there might be other restrictions on the final data that generate uniqueness. Note 
that the condition (21) is similar to the condition imposed in the one- dimensional case in Section 3. 
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